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This article reports results from numerical simulations of the gravitational radi- 
ation emitted from nonrotating relativistic stars as a result of the axisymmetric 
accretion of layers of perfect fluid matter, shaped in the form of quadrupolar 
shells. We adopt a hybrid procedure where we evolve numerically the polar 
nonspherical perturbations equations of the star coupled to a fully nonlinear 
hydrodynamics code that calculates the motion of the accreting matter. Self- 
gravity of the accreting fluid as well as radiation reaction effects are neglected. 



1 Introduction 

In this article we shall report of the work made in collaboration with J. A. Pons 
and J. A. Font (see Ref. P for further details). The purpose is to analyze the gravi- 
tational wave emission pattern from neutron stars as a result of accretion of matter 
which is shaped in the form of quadrupolar shells. Accretion is expected to happen 
following the gravitational collapse of the core of a massive star, once a neutron star 
has already been formed. Part of the remaining stellar material which has not been 
expelled by the shock driving the supernova explosion may fall back onto the neutron 
star, until a critical mass is exceeded and the star collapses to a black hole. Some 
more material may in turn form a long-lived, centrifugally-supported torus or disk 
if the collapsing star had initially some amount of rotation j2j . A detailed realistic 
modelization of the gravitational emission from accretion flows would require three- 
dimensional (magneto-) hydrodynamical simulations in general relativity, coupled to 
radiation transport and diffusive processes. However, some preliminary steps can be 
taken to understand the underlying basic physics in a qualitative way. The procedure 
adopted here lies in the borderline of full numerical relativity and perturbation theory 
and extends to relativistic stars the 2D-axisymmetric computations concerned with 
black holes, formerly performed in Ref. 0. The approach is as follows: the accret- 
ing matter, whose self-gravity is neglected, is evolved in a curved static background 
by solving the nonlinear hydrodynamics equations; the response of the star to the 
infalling matter, which triggers the emission of gravitational radiation, is computed 
by numerically solving the metric perturbations equations, with sources accounting 
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of the fluid motions. The key assumption of this method is that the mass of the 
accreting fluid is much smaUer than that of the star, so that, fluid self-gravity as weU 
as radiation reaction effects are neglected. The first approximation (no self-gravity) 
is in general valid for fluid motions in the vicinity of the star, where tidal forces dom- 
inate over the fluid self-gravity. The second approximation (no radiation reaction) is 
valid as long as the energy in the form of gravitational radiation is much smaller than 
the kinetic or internal energy of the fluid. We use then a particular implementation 
of the theory of metric perturbations of stars (see Ref. |4j and references therein for 
an extended description). A further simplification adopted is that the star is non 
rotating, so that metric, pressure, energy density and mass are obtained as solution 
of the Tolman-Oppenheimer-Volkoff equations of stellar equilibrium We assume 
barotropic equation of state written in polytropic form p = Ke^ to model a cold 
neutron star. We show results for two stellar models of mass M = IAMq with 
r = 2, the one being more compact than the other. The more compact model (A), 
has central density Cc = 2.455 x 10^^ g/cm^, and K = 122.25 km^, with R = 9.80 
km. The less compact model (B) has Cc — 0.92 x lO^^g/cm'^, and K = 180 km^ which 
leads to i? = 13.44 km. 



2 Even-parity perturbations of a relativistic star 

The presence of a fiuid distribution outside the star creates a small perturbation 
in its gravitational field, so that the total metric is expressed as g^i^ — g^^, + Sg^i,, 
where is the metric of the star and Sg^v is the perturbation; this can be decom- 
posed in odd-parity (axial) and even-parity (polar) modes, each carrying spherical- 
harmonic indexes I and m 0] . Since the hydrodynamics evolution is 2D-axisymmetric 
(9<p — 0) and because the star and the shell have no angular momentum, only the 
polar modes are excited by the accreting matter. For this reason we shall restrict 
the discussion below just to polar perturbations of stars. Formulations of the polar 
perturbations equations of a nonrotating stars suitable for time evolution are cur- 
rently available in literature (see Ref. pp and references therein). Our calculations 
are performed in the Regge- Wheeler gauge |H|, by specializing the gauge- invariant 
and coordinate independent formalism of Ref. . The equations obtained are equiv- 
alent to those of Ref. [H] and [J], although different metric variables are used. Thus, 
for each {l,m) pair, 6g^i, is parametrized by two scalar quantities, k (the perturbed 
3-conformal factor) and x (the actual gravitational wave degree of freedom), and it 
reads 

Sgf.^dxf'dx'' = [ix + k)e^''dt^ - 2V'e"+^didr + (x + k)e^''dr^ + kr^dil^] F^" . (1) 

On a static background ip does not play any dynamical role, as it can be obtained 
by quadrature from k and x l2j- We note that, despite this problem has actually two 
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Figure 1; Time evolution of the enthalpy perturbation of model B 



degrees of freedom inside the star and one outside it, we have found more convenient 
for numerical reasons to consider an additional variable inside the star, the perturba- 
tion of the relativistic enthalpy H = Sp/{p + e). We are left, then, with three degrees 
of freedom inside the star (x, k and H) and two outside (x and k). Outside the star, 
from k and x we get the Zerilli-Moncrief function ITTj 
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with X = l{l + 1), whose Fourier transform gives the energy radiated in GWs 
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Unlike other studies [51 E] , we choose to evolve the perturbations through a cou- 
ple of hyperbolic equations for x H plus an elliptic equation, the Hamiltonian 
constraint, which is solved for k at every temporal slice. With this "constrained" 
scheme, that works as a dynamical corrector of the truncation errors, we gain two 
important benefits: on one hand we can obtain the frequencies of the fluid modes 
(by Fourier transforming the waveforms) with very good accuracy even with a rel- 
atively low resolution; on the other hand, we do not have limits to the evolution 
time (unlike free hyperbolic formulations, where sooner or later a violation of the 
constraints occurs and the evolution must be stopped), and we can see dynamically 
the damping of the fluid modes due to GWs emission. Fig. ^ shows the evolution of 
an initial enthalpy profile a,s H = sin(7rr/2i?) for model B, observed at half of the 
star radius. The circles in the Fourier spectrum represent the fluid modes obtained 
using a frequency-domain code based on the Lindblom-Detweiler formulation of the 
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perturbed equations ^21 and described in The relative difference between the 
location of the peaks and these values is of some parts in a thousand, with a resolu- 
tion of 200 points. The right panel of Fig. prefers to the same evolution carried out 
for 2 s. The damping time of the /-mode for this model is t/ = 0.26 s: after ~ 1 s 
just the longer living p-modes survive in the signal ^ . 

3 Gravitational waves from fluid accretion 

We present an overview of a typical simulation of accretion onto model A. Ini- 
tially, the non spherical density profile is surrounded by a background fluid (an "at- 
mosphere" ) satisfying the stationary and spherically symmetric Michel solution . 
The initial rest mass density profile is given by 

P = PMich + Pmaxe^'^^''"''''^' sin^ , (4) 

with pMich being the profile consistent with the Michel solution. The mass of the 
shell is /i = O.OlAf , which corresponds to a maximum density of Pmax 3.5 x 
10-6 km^^ when we fix its width to k = 1. The (inhomogeneous) density profile of 
the atmosphere pMich is roughly three orders of magnitude lower than Pmax- The shell 
obeys a polytropic {p = ICp^) EoS with JC = 0.01 km^^"^ and 7 — 4/3. The internal 
energy profile is obtained from the baryonic density p and p through the first law of 
thermodynamics as e = p/ ((7 — l)p). The shell is initially at rest at tq — 20 km, and 
the waveform is extracted at robs = 250 km. 

The motion of the fluid in a curved spacetime is governed by the local conserva- 
tion laws of baryonic number and energy-momentum; as shown by |15| . the equations 
of general relativistic hydrodynamics can be written as a hyperbolic system of bal- 
ance laws with sources, and then solved using a Godunov-type scheme based on the 
characteristic information of the system. The hydrodynamics code employed here is 
the same of Ref. [3] . 

In a realistic scenario, such as e.g. fall-back accretion after a supernova explosion, 
the accreting matter interacts with the neutron star until a critical mass is reached 
and the star collapses to a black hole. As we cannot model this complex phenomenon 
within the current perturbative framework, we impose reflecting boundary conditions 
at the inner edge of the hydrodynamics domain (i.e. the surface of the star), so that it 
is seen by the external fluid as a hard sphere. This choice includes the most relevant 
effect, that is, the pressure gradient stops the infalling matter, and the accretion 
process is then followed by the formation of shock waves which propagate off the 
stellar surface. The impact of the shell perturbs the star and triggers its quasinormal 
modes of pulsation, so that it radiates gravitational waves. 

Fig. 121 exhibits the time evolution of the Zerilli-Moncricf function and the corre- 
sponding energy spectrum. In the waveforms, three phases can be identified. Firstly, 
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Figure 2: Waveform and energy spectra from a fluid shell falling on to star model A. 



the infalling, when the bulk of the shell is evolving outside the star, gradually ap- 
proaching it, which is characterized by the steady increase of the amplitude of the 
signal: it is very short as the shell is initially located very close to the stellar surface. 
Secondly, a burst-like peak appears in the GW signal, which, as found in simulations 
of gravitational core collapse , coincides with the moment when the shell reaches 
the surface, creating a shock wave which propagates off the surface. Finally, the ring- 
down phase, characterized by a GW signal which is not exactly monochromatic as a 
result of the complex interaction between the gravitational field of the star and the 
layers of fluid captured on top of the stellar surface in the process of readjusting them- 
selves to a new stationary solution. We note that, despite a remarkable resemblance 
with the results of core collapse simulations |EI , there are also important differences 
in the post-bounce phase dynamics, as the ringdown of the neutron star lasts for 
much longer times. In our idealized setup these pulsations are not quickly damped 
by the existence of a dense envelope surrounding the star, as happens in the core col- 
lapse situation. Concerning the energy spectra, the solid line is obtained by Fourier 
transforming the complete signal; correspondingly, the dashed line comes from the 
truncated waveform, in which we only take into account the contribution from the 
beginning of the burst (~ 0.9 ms) onward. The whole spectrum shows a complex 
structure: first of all, we find a narrow peak, identified with the /-mode of the star 
with an error around 1%. At higher frequencies, we have a broad band spectrum 
modulated by roughly equally spaced bumps {6f ~ 0.9 kHz). These are interference 
effects coming from the pre-bounce phase of the emission. This is confirmed by the 
fact that they are absent in the truncated spectrum. This latter presents one broad 
peak with a maximum at about 7 kHz, which is however too low to be identified 
with the w-mode. As mentioned before, its origin should be related to the motion of 
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the fluid shell and its interference with the gravitational fleld of the star: the high- 
frequency emission depends, then, on details of the dynamics of accretion rather than 
on intrinsic characteristics of the star. The total energy emitted in this process is 
computed by integrating the spectrum of Fig. |2| We get ~ 3.02 x 10~^ Mqc^ 
or, in terms of the shell mass, i?^" ~ 2.15 x 10~^ fic^. 

As a last remark, we note that a small amplitude peak, at frequency lower than 
that of the /-mode, is present in the spectrum. It is associated with oscillations of 
that part of the external fluid that has been gravitationally captured by the star as 
a result of accretion. We argue that the existence of this unphysical low frequency 
peak is an artifact produced by the reflecting boundary conditions imposed at the 
surface. In a realistic scenario, the accreted matter would not simply bounce at the 
stellar surface, but it would rather interact with the neutron star envelope, resulting 
in heating and suffering nuclear reactions, until it is reabsorbed by the star. 
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